June 1993 



ITP-SB-93-15 



Spectrum of Low-Lying Excitations in a 
supersymmetric extended hubbard model ^ 



Fabian H.L.Emer 

m 

and 



Vladimir E. Korepin 



Institute for Theoretical Physics 



^ State University of New York at Stony Brook 

^ Stony Brook, NY 11794-3840 

S ABSTRACT 
I 
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states. The spectra are found to be very different for repulsive and 
attractive on-site interaction. We also study the thermodynamics 
of the model. 
Pacs : 75.10.Jm 71. 20. Ad 



^This work was supported in part by NATO, S.P. Chaos, Order and Pattern (CRG 901098). 

^E-MAIL: FABMAN@MAX.PHYSICS.SUNYSB.EDU 
^E-MAIL: KOREPIN@MAX.PHYSICS.SUNYSB.EDU 



1. Introduction 

This is the fourth part of a series of papers deahng with the M(2|2)-symmetric 
model of superconductivity. In [1] a new model of strongly correlated electrons on 
a general J-dimensional lattice was proposed, which differs from the Hubbard model 
by moderate interactions of nearest neighbours. The additional interactions fall into 
two categories : the so-called "bond-charge repulsion" terms and nearest neighbour 
spin-spin and charge interactions are similar to the ones of the model of hole super- 
conductivity proposed in [2], and the t-J model respectively. The pair-hopping term 
(which allows local electron pairs to move coherently) has been studied extensively in 
the Penson-Kolb-Hubbard models^"^"^]. The model is of the "permutation type" first 
studied by Sutherlandt^l. The model has the interesting property that the ground 
states for attractive and moderately repulsive on-site interaction exhibit off-diagonal 
long-range order (ODLRO) and are thus superconducting in one, two and three di- 
mensions [''"^l. Furthermore it is possible to derive the exact expression for the unique 
ground state in the attractive regime in any number of dimensions. In one dimension 
the model is integrable and in [10] the nested Algebraic Bethe Ansatz solutions were 
constructed in detail. That paper also contains expressions for higher conservation 
laws and the proof of a lowest weight theorem of the Bethe Ansatz states with respect 
to the m(2|2) symmetry algebra of the model. 

This paper is a direct continuation of [10]. Starting from the Bethe Ansatz equations 
in the BBFF representation^ derived in [10] (see also [If]) we first derive the infinite 
set of coupled integral equations that drive the thermodynamics of the model (TBA 
equations). By means of taking the zero temperature limit of these equations, we de- 
termine the ground state of the model in one dimension. The result agrees perfectly 
with the result obtained in [10] by completely different methods. The zero temper- 
ature limit of the TBA equations also provides expressions for the dressed energies 
of the elementary excitations that are given by the Bethe Ansatz. We also study 
the limit of infinite temperature of the TBA equations. In section 4 we determine 
the spectrum of low lying excitations in the T = formalism. We carefully take into 
account that the Bethe Ansatz only provides lowest weight states of the m(2|2) sym- 
metry algebrat^°^. This makes it necessary to analyze not only the Bethe equations in 
the thermodynamic limit, but also to construct the complete m(2|2) multiplet from 
the lowest weight state obtained by means of the Bethe Ansatz. In order to enumerate 
all elementary excitations we use hoth Bethe Ansatz and global symmetry (see also 
[12-14]). A complete set of elementary excitations is of course necessary for deter- 
mining questions like charge and spin separation, or deriving effective field theories 

^The reason for choosing the BBFF representation is twofold : First it reduces to the Sutherland 
representation of the supersymmetric t-J model in the sector without localons. Second the ground 
state derived below is of a very simple form - there are no complex roots of the Bethe equations 
present in the ground state. 
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for the model. The results for repulsive on-site interaction (and filling greater than 
zero) fare as follows: Apart from spinons and holons, which are already present in 
the supersymmetric t-J model (which is a submodel of the M(2|2)-model) with linear 
dispersion relations around their repective Fermi surfaces, there exist electronic exci- 
tations with quadratic dispersions and a gap, and excitations involving local electron 
pairs. These also have a quadratic dispersion and a gap for // < (//is the chem- 
ical potential), but become gapless for // = 0. For attractive on-site interaction we 
find an infinite number of excitations with cosine-like dispersions. Amongst these are 
electrons and local pairs and their bound states. 

The hamiltonian on a lattice of L sites in the grand canonical ensemble is^^^ 

L ^ ^ L L 

H{tJ^, U,h) = H' + U Y,{n,,i - 7:){n,,-r --)-// J] n^- + 5](n^- 1 - n^- _i) (1.1) 
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where = — "^^-^ H^j+i is given by 



(1.2) 

Here U is the Hubbard model coupling constant, h is an external magnetic field, and 
/i is the chemical potential. By n- ^ = c\ ^c- ^ we denote the number operator for 
electrons with spin a on site i and we write n- = n--^-\- n- Under the particle-hole 
transformation c^^ ^ the hamiltonian H{ji^U^h) transforms (up to a constant) 
into H[—fi^U^h) so that it is sufficient to consider only the region // < of the 
chemical potential. It was shown in [1] that is invariant under a m(2|2) symmetry 
algebra. The three other terms in H are elements of the Cartan subalgebra of m(2|2) 
(and commute with and each other), and thus and H have a complete set of 
simultaneous eigenstates. The dynamics of the model are such that the individual 
numbers N.^ and of electrons with spin up and spin down, and the numbers Nj 
and Nf^ of doubly occupied ("local electron pairs") and empty sites ("holes") are 
conserved quantities. We will choose the following conventions throughout this paper 

N^- = number of single electrons with spin up 
= number of single electrons with spin down 
= N^- -\- = number of single electrons 

Ni = number of local electron pairs 
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Nf^ = number of holes 

~ ~^ — number of "bosons" . 

Let us now review a few facts about the u{2\2) symmetry of the modeP, which 
we will need in section 4 below. There are a total of 16 generators, 8 of which are 
fermionic and bosonic respectively. The eight bosonic operators fall into two su{2) 
and two uil) subalgebras. The two su{2) algebras are generated by the spin-operators 



i=i i=i i=i 

and by the ?y-pairing like operators^^^^ 

L L L ^ ^ 

i=i i=i i=i 

The two charges are given by the identity operator and the Hubbard interaction 

^ = EKi-^)K-i-^) • (1-5) 
The eight fermionic generators are given by 

L L 
Qa = - ^J,-a)Cj,a , Qa = J^^J.-aS,^ , (7 = ±1 (1.6) 

and their hermitean conjugates. The lowest weight theorem states that all states con- 
structed by Bethe Ansatz are lowest weight states of the u[2\2) symmetry algebra^^^l. 
In the BBFF representation of the Bethe Ansatz this implies that {\BA) denotes a 
Bethe state) 

r^lBA) = = S\BA) = Qi\BA) = QJBA) . 

A complete set of states can be obtained by acting with the u[2\2) raising opera- 
tors on the Bethe states. Every Bethe state gives rise to one u[2\2) multiplet. These 
lowest weight representations fall into two categories: so-called typical and atypical 
ones^^^'^'']. The dimensions of the typical representations are given by 



dim(iVpiV^,iV, 



= 16(iV^ + 1){L -N^-N^- 2Ni + 1) , (1.7) 

typ 



where — A^^ -|- 1 and L — N-^ — N ^ — 2Ni -\- 1 are the dimensions of the representaions 
of the spin and ?y-pairing su{2) algebras respectively, and the factor of 16 is due 



^More precisely of H'^ , as only H'^ commutes with all generators of the m(2|2) algebr 
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to all possible actions of the fermionic generators. For the atypical representations 
some (nontrivial) actions of the fermionic generators give zero, and the multiplets are 
smaller. Two such cases appear for the model at hand 

dim(iv,0,0) =4iv, 

(1.8) 

dim(iVp 1,0) = dim(iV^ - 1, 0, 1) = 8{N^ + Ni-1){L-N^- iV;) + iL . 

Whereas (1.2) commutes with all 16 generators of m(2|2), the hamiltonian (1.1) does 
not as the chemical potential, magnetic field, and the Hubbard interaction break the 
symmetry in a trivial way. The energies of the m(2|2) descendants of the Bethe states 
\BA) are shifted as follows (for zero magnetic field) 

H{fi,U,0) \BA) = E \BA) 
i7(/i,[/,0) \BA) = E \BA) 
i7(/i, [/, 0) r^^ \BA) = {E- 2/.) r^^ \BA) 

H{^i, U, 0) Ql \BA) = {E + ^-iJ^)Ql \BA) 

H{iJ^, U, 0) \BA) = {E + ^ + iJ^)Q^ \BA) . 

Equations (1.9) will be important for the classification of the excitation spectrum in 
section 4. 



2. String Hypothesis and Logarithmic Form of the Bethe Equations 

In this section we classify all roots of the Bethe equations in the BBFF repre- 
sentation[^°'^^] by means of a usual "string picture". After inserting this classification 
into the Bethe equations and taking the logarithm, solutions of the Bethe equations 
are unambiguously parametrized by a set of integer of half-odd integer numbers. This 
is the starting point for both the Thermodynamic Bethe Ansatz constructed in sec- 
tion 3, and the determination of the excitation spectrum at zero temperature (section 
4). 

The Bethe Ansatz equations for the BBFF grading are of the form 

_TJ ^k-^ +i i-r A,-Afc + 2z 
K + ^) ~l\K->^^-^ l\ A,-A,-2z ^ " 1' " " " ' + 



k=l ^-j « „ = i ^ra ^.j I 

n \ -AL^ + 2» Y[ -1 - rn-\ N 

A«-A«-2z 11a«-A;.^) + z ' 

1=1 ' J=l J 



f2.11 
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Here we have rescaled the spectral parameters by a factor of 2 as compared to formula 
(2.78) in [10]. Energy and momentum of a state parametrized by a solution of (2.1) 
are given by (in the grand canonical ensemble) 

KA„...,AJ = (Z + l)7r-z J2 HyT^) ■ 

j=i J 

(2.2) 

In the limit of large lattice lenghts L the imaginary parts of solutions of (2.1) 
are distributed according to the following idealized presciption (string hypothesis) 

A-strings : A^-'' = A^ + z(n + 1 - 2j) , A^ G M , j = l...n , (2.3) 
real A^'^ : Af;) G M , (2.4) 

r AW^^' = A«+z(3-2j) , J = 1,2 
A(^)-A(^) -strings , (2.5) 
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A(^) -strings : Xi'^"" = X^'^" + t{p + 1 - 2k) , X^'^" G M, k = l...p . (2.6) 

Inserting (2.3)-(2.6) into (2.1) and "multiplying out the strings" we arrive at the 
following equations involving only the real centers of the strings 



m=\ n n n 

\(f3,m):^(a,n) J 7=1 5=1 

n n <^)= n n k^^m^) n ka^^^^-^/^^) 

(2.7) 

n n^-(^^/^^-^^)= n <>^''->^n n k^^^) 

ri=l 0- = ! 7=1 e=l 

7=1 (a',05^(''>p) 

Here e(x) = ^ and 
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Y ^ ) - )...e^{ - )ef^^) if n ^ m 

V |ri — m| / V |ri — m| + 2 ' ' ' ' V n-\.^^ — 2 ' V n-\-m ' ' 

K.M = { (2.8) 

are Takahashi's e-functionst^^'^^l, and 

(i) denotes the number of A-strings of lenght n 

(ii) M^^"^ denotes the number of real A'^^^'s 

(iii) Af2^^ denotes the number of A*^^^ — A'^^^-strings 

(2) 

(iv) Mp denotes the number of A'^^^-strings of lenght p. 
Furthermore we define for later convenience 



M:=Y,nM^ = N^ + N, + N, , 

n=l 

M(^) := M^^^ + 2M^^^ =Nj^ + N^ , 

CO 

■p=i 

Taking the logarithm of equations (2.7) we arrive at 



(P;in)^(a,n) 7=1 5=1 

E E ^(^)=2vr V 5: 5: ^(^^M^)+ E ^(A^^)^-A(/)^) 

n=l a = l 7=1 e=l 

^^(&^)=27riV,^+ E ^,,(aP)^-aP)') . 

7=1 (ii,imi',v) 



(2.9) 



Here J^, K^, are integer or half odd integer numbers, 6{x) = 2arctan(a;), and 
^(^S) + 2^( i "n, ) + --- + 2 g(— ^) + g(^) ifn^^m 

V |ri — m| / V |ri — m| + 2 / V ri+m — 2 / V ri+m / ' 



2^(f) + 2^(f) + ... + 2^(^) + ^(|^) ifn = m . 

(2.10) 
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The integers (or half-odd integers) are determined by the choice of cuts in the 
logarithms and have the following ranges 



J. 



integer 

half — odd int. 



if < 



CO CO 



Ml - m\^^ - M^^^ - 1 



even 
odd 



f2.ii: 



CO 

m = l 

CO CO 

n=l p=l 

CO 

\Ks\<-iM, + 2Y,M^-M['^-Mi'^-l) , 

m = 2 

CO 



(2.12) 



q = l 



Here = 2 min{m,n} — S^^. These equations will be needed to determine the 
ground state and the excitation spectrum at zero temperature. 

Insertion of (2.3)-(2.6) into (2.2) leads to the foUwoing expressions for energy 
and momentum in terms of the (real) centers of the strings 



TT oo M„ , 

n=l a = l ^ 



4n 



+ 2hn 



U 



+ ikff V +l7 + h) + M^^\U + 2h)-Y, 2/«P Ml 



2) 



f2.13l 



n=l a = l 



n 



Inserting the logarithmic form of the Bethe equations (2.9) into our expression 
for the momentum (2.13), it is possible to express the momentum completely in terms 
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of the integers present in (2.9) 



oo M„ 



, n 

n=l a = l 

f ( E E + E + E - E Z I + ' ( i + 1 - E 

ri=l 0- = ! 7=1 6=1 p=l 1^=1 j ri=l 



(2.14) 



3. Thermodynamic Bethe Ansatz 

In this section we will derive an infinite set of coupled integral equations (TBA 
equations) that determine the thermodynamical properties of the model in terms of 
the values of the magnetic field, chemical potential and coupling \J . We will follow 
very closely the route taken by Takahashi in [18, IQ]"*, which is a generalization of 
Yang and Yang's methodt^^^ to the case of an infinte number of string solutions of the 
Bethe equations and the nested Bethe Ansatz. Before we start with the derivations 
we recall that it is sufficient to consider values of the chemical potential < 0. One 
might expect problems with the formalism at zero chemical potential for the following 
reason: It was shown in [7] that for // = the compressibility becomes infinite, and 
the ground state structure is governed by the m(2|2) symmetry, which is a structure 
that complements the Bethe Ansatz rather than being encompassed by it. However, 
our analysis at zero temperature shows, that taking the limit // — )• of the TBA 
equations derived below (from the region // < 0) gives the correct results. This leads 
us to believe that the TBA equations also give the correct thermodynamics at // = 
for r > 0. To derive the TBA equations we first rewrite the Bethe Ansatz equations 
(2.9) using the concept of counting functionst^^^. Defining z,^{\) = 27r-^ (and similarly 
y^{\^^'') = 27r-^ ,etc.) and then taking the thermodynamic limit Z — )• oo in which we 
replace sums by integrals weighted by the density functions we arrive at 



^n{^)= ^ ii) - >^ / e.^.^{X-A)p.jA) 
-1 



+ 00 /' + 00 



dAe{^)p^^\A)- I dAe^^,{A-X)p^i\A) 

CO 

j/,(aW)=5] / dA9i^)p^^ 

_-|^ — oo 

o 

Oi^)pi^\A) - I .A ^(a(^)-a)^^,^^(a) 



(3.1) 



OO 



p=l 



+ 00 /' + 00 



oo 



^For a review see [20]. 
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1 'J —oo 



n=l 



+ 00 p-\-oo 

.A ^(a(^)-a)^(^)(a)- / .A e{^^)pf{K) , (3.2) 

■00 J —CO 

/ + 00 /' + 00 

.A ^(^)^(^)(a)- 5] / .A ^,,(A«-A)/>r (a) • 
■00 q = l ^ 

Note that in the above equations all spectral parameters are allowed to go to 
infinity, which corresponds to the action of the symmetry algebra m(2|2). Thus the 
Thermodynamic Bethe Ansatz takes into account not only Bethe states (which are 
only the lowest weight states of m(2|2) multiplets) but also states obtained from those 
by acting with m(2|2) raising operators. This is absolutely essential because all states 
in the Hilbert space contribute to the thermodynamic properties and have to be taken 
into account. 

The densities for the vacancies (allowed integer numbers parametrising solutions 
of the Bethe equations in the sense of (2.9)) are given as usual as the derivaties of 
the counting functions with respect to the spectral parameters 



(3.3) 



2vr(^(,^^(A(^)) + ^r(A(^))) , 
2vr(^(,^^(A(^)) + ^t^(A(^))) , 

2vr(,(^)(A(^)) + ,tHA^^^)) • 

Here the /j's are the distribution functions of spectral parameters corresponding 
to the taken integers ("particles") whereas the /j's are densities related to the allowed 
integers which are not occupied ("holes"). Defining the integral operators 



JA(i) 

JA(i) 
Jxp(A(^) 



N */ 



1 d 

= ^n-r.fW + ^^ JA'0„,„(A-A')/(A') , 

^ (3.4) 



and differentiating (3.f),(3.2) with respect to the spectral parameters the Bethe 
equations (in the thermodynamic limit) turn into 
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1 2n 



27r n2 + A2 



n.m r m 



m = l 

CO 



+ (^.,l-'^.,l)*/'2 



(1) 



,(1) 



oo 



n=l 



? 

OO 



A 






OO 


A 




A 


-[1]*P'2' 

A(i) 


A(i) 


-E 




A(i) 






A(i) 


-[1] 


*p['' 


- [2] * P'i' 

A(i) 


A(i) 



,(2) 
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p=l 



+ [P] * /O^l'^ 



A(2) 



A(2) 



(3.5) 

This set of equations expresses the densities for the "holes" in terms of the 
densities for the "particles". The set of equations determining the densities p are 
obtained from the condition that the system be in thermodynamic equilibrium, i.e. 
from minimization of the thermodynamic potential 

F 1 U 

OO /i i-yQ 

-J2 ^^{(1^ + 2-'') Pn + T{{pn+Pn)HPn + Pn)-PnHPn)-PnHPn))} 
n=l •' 

/OO 
jA(^) {(.+f +.) - T iif^^ + Hpi' + Pi') - Pi' Hpi') - p['' Hpi')) } 
-00 

+ / JA(^) {iu+2H) fJi^ - T {{p[^^ + ^tO Hp'2 ' + pi') - pi' Hpi') - pi' Hpi')) ] 



■p=i 



J2 I dX^'" Wp p'p' + T {{l^P + 1^^) + 1^') - P'f Hp'p') - P'p' Hp'p')) } • 

) 

(3.6) 



At this point it is convenient to introduce the quantities 

P 



a 



(3.7) 



_ - fi - Pli_ 

1 (1) 1 P2 (1) 1 Ip (2~ 

Pn Pi P2 PP' 

Varying the potential F under the constraints (3.5) we obtain the following 
conditions for thermodynamic equilibrium 



1/^4 \ 00 y ^ 



N*ln( 1 + ^ 



-iA,i-K,i)*'i^ (1 + -^ 

A V P2 



(3.8) 
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ln(/3,(A(^))) 



T 



+ [l]*ln( 1 + -^ 



/I 
*ln 1 + — 

ri=l ^ ' 

- Vb]*ln 1 + - 

1) V 7, 



A(i) 



(3.9) 



ln(/3,(A(^))) 



U + 2h 
T 



oo . 
n=l ^ 



A(i) 



+ [2]*ln(l + i- 



+ [l]*ln (1 + -^ 

A(i) V Pi 



A(i) 



(3.10) 



ln(l+7,(A«)) 



0=1 ^ '1 



+ b]*ln(l + -^ 

A(2) V Pi 



The thermodynamic potential at equihbrium takes the form 



F 



U 



L ^+^-^ 4 



00 

rV- 

^ 27r 



In 1 + 



2n 



A(2) 



(A) y n2 + A2 



(3.11) 



(3.12) 



The TBA equations (3.8)-(3.12) determine all thermodynamic properties at finite 
temperature as functions of the chemical potential the magnetic field h and the 
coupling U. 

Using methods similar to the ones employed in [22] it is possible to recast (3.8) 
and (3.11) in "tridiagonal form" 



ln(a„(A)) =p*ln (l + a„_i) 
/ 47r 



+ P*ln(l + a„+i^ 
1 



\4:T cosh{^' 
ln(7,(A«))=p*ln (I + 7.-1) 



+ p * In 1 + 



/3i 



•^n,! 1 + 



/32 



-'n,2 1 



A(2) 



+ ]5*ln (1 + 7,+!^ 



+ ]5*ln(l + -^ 

A(2) V Pi 



Ss,l- 



A(2) 



(3.13) 



Here p is an integral operator defined by 



P* f 



A J -00 



dX'- 



4cosh 



((A-A')f) ^^^'^ 
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We have not succeeded in simplifying (3.9) and (3.10) in a similar manner, which 
makes it nontrivial to take the T — )• limit of the TBA equations (see below). Similar 
problems have also been encountered for other models^^"^]. Equations (3.8)-(3.11) can 
be used to compute the low-temperature specific heat^^'*]. 

The limit T^O 

In order to study the zero temperature limit we first define the dressed energies 
according to 



e e[''^ 

«n = exp(^), /3i = exp(^), (3^ = exp(^), 7^ = exp(^) . (3.14) 

We will show in our discussion of the excitation spectrum in the T = formalism 
below that the e's really play the role of dressed energies at zero temperature. It can 
be shown that this interpretation also holds at finite temperaturet^^'^^l 

From (3.13) it is clear that all e^/^ and all e^, n > 3 are non-negative and will 
drop out of the r.h.s. in the T — )• limit of the TBA equations (3.8)-(3.11). The only 
energies that can change sign are e^^, £3, e'--^^ and £2^''. Investigation of the TBA equa- 
tions shows that either both £2^'' and £2 change sign, or both of them are non-negative 
for all values of spectral parameters. By a completely independent computation in 
the r = formalism we can show that for T = both are nonnegative for all values 
of spectral parameters. In the light of this we conjecture that even for small positive 
temperature only and e'--^^ change signs as functions of their respective spectral 
parameters. The ground state in the zero temperature limit is determined by using 
the knowledge which of the dressed energies are non-negative to simplify (3.5). We 
find that in the limit T — )• only the densities p^iX) and p'--^^ {X'^^^) are nonvanishing. 
It can be seen that Pi{X) is positive for values of A in the interval [—A, A] ( where 
A is some positive real number), and vanishes outside this interval. A similar result 
holds for p'--^\X'^^^) and [—B^B]. The values of A and B depend on the values of the 
chemical potential the magnetic field h and the coupling U. In the limit of van- 
ishing external magnetic field (/i — )• 0) A tends to infinity, which is the case we will 
restrict ourselves to for most of the remainder of this paper. Using (3.5) we thus find 
that the densities describing the ground state for p < are given as the solutions of 
the following coupled integral equations 

12 I 4 

^^^^^ = 2^TT^ - 2^ y_ ^^'4 + (A-A')^ P^^^'^ 

+ — r^'jA^^) ^ , p''\X'''^) , (3.15) 

27r J_B 1 + (A- A(i))2 ^1 ^ ^ ' > 

1 1"^ 2 
^W(A(^)) = — / dX p, (A) . 

^ ' 2^ J -00 1 + (A-AW)2 ^ 
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The value ol B is determined by the requirement that 



N N 

JA(^Vi'^(A^'^) = -^ = l--^ = l-D. (3.16) 
-B L L 

In the grand canonical ensemble B is given as a function of the chemical potential 
and the coupling U via the boundary condition e'--^^{±B) = 0. This enables us to solve 
for D as a function of 

2 ln(2) 



Q 




.5 1 1.5 2 

/J, + U/2 

Figure 1: T = Ground state density as a function of cliemical potential and 

coupling U for < 

Equations (3.15) are recognized to be identical to the ground state integral equa- 
tions of the supersymmetric t-J model (in zero external magnetic field) in the Suther- 
land representationt^'^^'^*']. An investigation of the dependence of the integration 
boundary B on the chemical potential fi and coupling U shows that the ground state 
of the u[2\2) model for // < is actually identical to the ground state of the super- 
symmetric t-J model for a value of the t-J chemical potential /i^j = // -|- ^. This 
confirms a result stated in [7] for the u[2\2) model in any dimension. From Figure 1 
we see that the ground state is equal to the empty vacuum |0) in the region where 
/i -|- ^ < 0. For /« -|- ^ > 21n(2) the ground state is equal to the half-filled ground state 
of the supersymmetric t-J model. Thus we exactly reproduce the // < part of the 
zero temperature phase diagram derived in [10], which is pictured in Figure 2. 
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st-J half filled band (insulator) 




IV ^--^ 


III |st-J) 






*- II ^^^^ 


— empty state 


filled state — 


V 




^ 1 





Figure 2 : Ground states in the grand canonical ensemble 



Taking now the T — )• limit in (3.8)-(3.11) and using (3.14) we arrive at the 
following equations for the dressed energies at zero temperature and h = 

1 2n 



CO 



2 27ri_^ i + (A-AW)2 

' ^ ^ 27ri_^ 4 + (A-AW)2 1^ ^ 

1 /'+^ 2 



(3.17) 



2. ^^TTp^^ ^"(^^ ' 

-^^^(^")--^^^-2i/r-^ v+(A''-A)- -"^^^ 

In section 4 below we will reobtain these equations as describing the excitation 
energies of the low-lying excitations above the ground state in the zero temperature 
formalism. The integral equations for the dressed energies cannot be solved analyt- 
ically for arbitrary filling. However one can decouple them by Fourier transforming 
and finds 
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27r I r . 2 1 



u.. /-^ 1 



( A) = ^ + _ _ 47r i?( A) + / JA' i?( A - A') e^^^ (A') 
where = ^ /_";;^ i+(J_,)2 ^^^iW. 



The limit T ^ oo 

We now turn to the study of the infinite temperature hmit. The analysis is 
a straightforward generalization of Takahashi's method for the Fermi gas with 6- 
function interaction^^*]. In this section we will keep the possibility of having a finite 
magnetic field h. The infinite temperature limit is of interest because it gives an 
estimate for the number of Bethe states in the thermodynamic limit, thus providing 
a rough check on whether the Bethe states form a basis of the complete Hilbert space. 
Furthermore it is possible to obtain a high temperature expansion for the Gibbs free 
energy in this limit. 

We first observe that for T — )• oo all driving terms tend to zero in the integral 
equations and a^, /32, 7^ all turn into constants. Therefore the integral equations 
turn into coupled algebraic equations, which can be solved analytically. We first note 
that in the infinite temperature limit insertion of (3.8) for n = 1 into (3.10) leads to 
the following simple relation , which allows us to eliminate from all other equations, 
which we will do below. 

/3 = — ^ — 
The infinite system of coupled integral equations (3.13) turns into 



'1 + - 

^ /?! 



(7i)^ = (i+72)(i+-^: 



+ (7j' = (l + %-i)(l + %+i), s>2 



(«n)^ = (l + «n-i)(l + «n+i), ri>3. 



(3.19) 
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and we still have to consider the equations for a^^ and 



Mi + «.) = 2gl„(i + i-)-l„(i + l 

Using (3.19) in (3.20) we find the "boundary conditions" 



fi-\- h 
T 



(3.20) 



lim In = 7^ , In (/?! ) = - - In ( 1 + 7i ) 

limln' ^^ + ' ^-'^ 



(3.21) 



The recursion relations (3.19) have the solutions 



7. = {fis)r - 1 , /3i ^ 



(/(o))^-i 



(3.22) 



where 



= 

Conditions (3.21) now determine a, 6, and c as functions of T, // and /i (Note 
that like Takahashi we keep the ratios ^ and ^ finite for the time being, but take ^: 
to zero.) 



cosh 



(/32)-' 



cos/ 

\ 1 / 

2 



/ szn/i (2^ , , 
a = - 1 , n > 2 



/3i 



1 



(3.23) 



4:{cosh{^)) -1 



sink (^) ' 
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The thermodynamic potential F can be evaluated exactly in the infinite temper- 
ature limit. Clearly the entropy contribution is dominating and the result thus gives 
us an estimate for the total number of states in the Hilbert space. This is because 



- = -- In (tr{e^^^^)) > - - ln{tr{id)) = -rin(4) . 

Li h \ ' T^oo J-j 

The formalism developed above leads to the following result 



| = l + /,_^_|_r^ln(^l + ^^^ -Tin (^2 cosh (^) + 2 cosh (^^^^ . 

This obviously approaches — rin(4) and we thus conclude that the dimension of 
the Hilbert space is roughly 4^ (because there can be corrections of order In(iv) which 
we cannot determine using the techniques employed above), which is the correct 
result. A rigorous clarification of the question whether the set of Bethe states is 
complete is rather difficult, as counting arguments as employed in [14] (for the XXX 
model) or [13] (for the Hubbard model) fail due to the presence of zero-norm statest^^^^. 
Thus it is necessary to first determine the norms of the wavefunctions["^°] and then 
count the non- vanishing statest"^^!. 



4. Excitation spectrum at zero temperature 



In this section we will derive the spectrum of the low-lying excitations over the 
ground states in the various regions in the fi — U plane at T = 0. The external mag- 
netic field h is taken to be zero. We will make use of the "standard" T = 0-formalism 
that goes back to Hulthen and des Cloizeaux and Pearsont"^^'"^"^!, the results of which 
are found to agree with the T — )• limit of the thermodynamics as discussed in section 
3. In [10] we have proved, that the Bethe Ansatz only provides lowest weight states 
of the m(2|2) symmetry algebra. Therefore the spectrum of the M(2|2)-symmetric part 
of the hamiltonian (1.1) is highly degenerate. The chemical potential and Hub- 
bard interaction terms break most of this invariances. The only degeneracies left are 
due to the spin-SU(2), which is not broken as the magnetic field is taken to be zero. 
In order to construct the low-lying excitations over the ground states (which were 
determined in [7], and again from the TBA in section 3 above) we will first use the 
usual Bethe Ansatz techniques to obtain the excitations given by the Bethe Ansatz. 
We will then construct the complete m(2|2) multiplet of elementary excitations by 
acting with the m(2|2) lowering operators. 



This means that the wave functions to certain solutions of the Bethe equations vanish. 
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4.1. Ground State and Excitations in Region III 

Region III in Figure 2 is the region where /i<00</i + f< 21n(2). We nave 
shown is section 3 that the ground state in this region is given by the ground state of 
the supersymmetric t-J model for the induced value of chemical potential /i^j = // + ^. 
We will now proceed to construct a complete set of excitations over this ground state. 
Naturally we will encounter all excitations that are already present in the supersym- 
metric t-J modeP. In addition to these we will find a large number of "localonic" 
excitations, that is excitations involving local electron pairs. We will start with a 
discussion of the ground state and then give a complete account of all t-J-like exci- 
tations. The BBFF representation of the Bethe equations^^^^ we have chosen in this 
paper reduces to the Sutherland representation^^] of the supersymmetric t-J model 
in the sector without local electron pairs. The ground state and excitation spectrum 
of the t-J model in this representation was studied to by Bares, Blatter and Ogata 
in [26] (see also [34]). They studied the spin-singlet and spin-triplet excitations, and 
the elementary charge excitations. Forster and Karowski^^^l completed this analysis 
by taking into account also the m(1 |2)-descendants of the Bethe states. Below we will 
perform a similar analysis for the m(2|2) case. In order to discuss the m(2|2) structure 
we will need the commutation relations (1.9) between the hamiltonian (1.1) and the 
symmetry generators of the m(2|2) algebra. 

Ground state 

The derivation of the ground state and the excitation spectrum in the zero tem- 
perature formalism is based on equations (2.9)-(2.12). In the BBFF representation 
and the region // < the ground state is given by the Bethe Ansatz (In the BFFB 
representation, which reduces to the Lai repr. of the supersymmetric t-J model in 
the sector without local electron pairs, the ground state is not given by the Bethe 
Ansatz. It is obtained by acting with Q\Q^_i on the Bethe state in the multiplet. 
This action is equivalent to sending two spectral parameters in the Bethe state to 
infinity, which is how the ground state in the Lai repr. was constructed by Bares et 
al.). It is the lowest weight state of a m(2|2) multiplet of dimension 16. The ground 
state is in general degenerate with other Bethe Ansatz states^ but there are two cases 
where it is not : 

(i) L even, Nf^ even, and odd 

(ii) L odd, Nf^ odd, and odd. 

We will consider only case (i) for simplicity. 

The 15 other states in the multiplet are not degenerate with the ground state for 
less than half filling (as the supersymmetry is broken by the chemical potential and 

^This is because all eigenstates of the supersymmetric t-J model are also eigenstates of the m(2|2) 
modelt^]. 
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Hubbard interaction according to (1.9)), so that there is a unique ground state. In 
terms of the distributions of roots of the Bethe equations it is described by taking 



Ml = iVf ^ + iVf ^ , M„ = Vn > 2 , M^^^ = iVf ^ , M^^'^ = = M^^) , 

and by completely filling all vacancies for between — /^^^^ = ^{Nh '^ + ^f'^ ~ ^) 
and /^(j^ (symmetrically around zero), and by filling all vacancies for symmetrically 

1 i T iV«^-l 7 

around zero up to J^^^ = —'^ — . 

Inserting this prescription into (2.9), subtracting equations for consecutive in- 
tegers, and then taking the thermodynamic limit, we arrive at the set of coupled 
integral equations (3.15). 

The ground state energy is computed from (2.13) (taking into account that /i = 0) 



-.■.l + h-f^-^+ r dX gi(A) pf'iX) + dA g^(A) p['''''{X) - 



(4.1) 

where ^|(A) = —^2^ and ^2(^) ~ + ^ "bare energies" of the ground 

state particles in the A- and A*^^^- distributions respectively. In the second line we 
consider the ground state energy as a function of the integration boundary i?Q, which 
will be convenient for the analysis below. 



Spin-Wave Spectrum 

The spin- wave-type excitations are very similar to spin- waves in the spin-| 
Heisenberg XXX antiferromagnet. The correct elementary excitations for that case 
were identified by L.D. Faddeev and L. Takhtajan in [14,35]. We will follow their 
analysis very closely. We start by trying to insert one "hole" into the distribution of 
integers i.e. choose M^ = iVf ^ + Nf^ - 1. The number M^^^ of A^'^'s is kept fixed 
at N^'^ . This however changes the number of vacancies for to N^'^ -\- N^'^ -|- 1, so 
that there are really two holes. We denote the spectral parameters corresponding to 
the missing integers by A^ and A2. By a similar analysis we find that we only can 
make an even number of holes by this procedure. In the thermodynamic limit we find 
the following set of coupled integral equations for a state with 2N holes at positions 
A^ j = 1...2iV 

In what follows we denote the number of holes and down spins in the ground state by #f and 
NP^ respectively. 
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+ 2^ "^'"'tTV^ .n.V)-^E^(A-Aj) , ,4.2) 



1 f°° 9 
«W(A(^)) = — / d\ , pJX) . 



In the next step we now introduce a matrix integral equation formalism which 
will be very convenient in deriving the excitation energies. We first define the integral 
operator K 



(4.3) 



A^ii(A) = - K,,{X) = = K,,{X), K,,{X) = . 

Using the vector notation 

9. = {p^,P^^r, '^l=(lTTX^^^y (4-4) 



we rewrite (4.2) as 



^a, - ^K,) Qp = Z + ^ , (4.5) 



where the "order 1 corrections" t/j^ are given by 



The excitation energy E^^ = E[B) — Eq^^Bq) follows from (2.13) to be 

-B 

The differences of densities 



E.. = L j d\f,{X) {p,{X)-pf'{X)) + L I JAg^(A) (^(^)(A)-^(^)^"(A)) . 

(4.6) 



(4.7) 

a,(A) = Z(^(^)(A)-^(^)^"(A)) 
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obey a similar set ol coupled integral equations as the densities themselves 



^afj - ^K^f, ]af, = i^^ . (4.8) 



The inverse of this equation reads (assuming invertibility) 

<Jp=Ud-^K] V„- (4.9) 
The dressed energies at zero temperature are defined according to 



<^„/3-iO^/3 = ^l' (4-10) 



27r 

where are defined above and we impose the boundary conditions (vanishing 
on the dressed energies at the Fermi boundaries) 

e^(±oo) = 0, e2(±5) = 0. (4.11) 

This definition of course coincides with the integral equations for the dressed 
energies derived from the thermodynamics if we identify ei(A) = e;^(A) and e2(^) — 
6^\tC). We now can rewrite the expression for the excitation energy using (4.3), (4.6), 
(4.8), (4.9) and the inverse of (4.10) 

2 2N 

where we have used the expressions for t/j^ in the last equality. A numerical solu- 
tion of the coupled integral equations for e-^ and e'--^^ shows that both dressed energies 
are gapless at T = 0. Therefore the (multi-parametric) excitations constructed above 
are gapless as well. For a plot of the resulting dispersion we refer to [26]. We see that 
there is no charge associated with this type of excitation, whereas the component of 
the spin in the quantization direction and the total spin are seen to be 



S' = N , S = N. 

The value of the total spin follows from the value of and the lowest weight 
theorem. The simplest example is the spin triplet excitation constructed in [26]. The 
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momentum of the excited states can be evaluated by means of (2.14). For even N it 
is simply given by 



2„ 2JV 

P = V /'^ mod 27r . 

L ^ ^ 

1=1 

For odd N the integers cannot be distributed symmetrically around zero, and 
there are two branches of the excitation (depending on whether we fill the "integer" 
— ^ or — ;r— ). The momentum is 



P = ^(Z±iVr)-f 5:^t = ±2^^-f 2., 

where the Fermi momentum is defined as kp = 

It is clear that the type of excitation contracted above does not exhaust all 
possible spin-wave-type excitations, as for example the spin-singlets are not included. 
In order to construct these we consider a state with 2N — 2 holes in the integers 
and one A-string of length n (with integer /]^), i.e. 

M, = N^' + N^'-N , M„ = l, M„ = Vm^{l,n}, 

^{1) ^ ^GS ^ ^(1) ^ Q ^ ^(2) _ 

As compared to the "pure" hole excitation above we have an additional n-string, 
which "eats up" two vacancies in the /^-distribution, so that we are left with 2N — 2 
holes. The integer for the n-string can be computed using (2.12) and one finds that 
l-^a I ^ {N — n) are the only allowed values. The spectral parameter describing the 
center of the n-string can be determined from (2.9) 



LO {- 



L 



/3=1 



M 



(1) 



=1 



0{ 



.(1)1 



In principle this equation fixes A^ in terms of the A^'s, but the precise connection 
is not important for our purposes here®. 

Taking the thermodynamic limit of the Bethe equations (2.9) for this type of 
excitation we obtain the following set of integral equations for the densities describing 
the distributions of rapidities 



For n = 2 and # = 1 it follows from a computation like in [26] that \\ = 
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12 1 ... 4 

2^1 + A2 ~ 2^ 



1 2 

27r 7_R 1+ A- AW 2 ^1 ^ ^ 



2JV-2 



lye,. ... 1 / 2(n-l) I 



2(n + r 



-)' + (A-A?)' (n + l)2 + (A- A?|')2 

(4.12) 

Defining 



9 A/" 

/ Sr,s(\ AM ^ 2(n-l) 2(n + l) 

■'■^ 27r V(?^-l)2 + (A-A-)2 ^ (n + l)2 + (A-A-)2^ 



1 2n 

V2 



2tt n^ + {\-Xlf' 

we arrive at a matrix equation for the backflow distribution functions ^ (defined 
in (4.7)) which is of the same form as (4.8). After similar manipulations as for the 
case with only holes we arrive at the following expression for the excitation energy 

2^ /""a 

E,, = gi(A^) + Y, d^a eM i^a , (4.13) 

where ^"^(A^) = — ^2_(_'^("n-|2 is the "bare" energy of the n-string, and are defined by 
(4.10). Recall that e-^ 2 is a matrix notation for e-^ and e'--^^ 

This expression breaks up into contributions from the holes and the n-string. 
Defining the dressed energy e^(A^) of the n-string according to 

+ 27r n2 + (A(^) - A)2 ^ ' 

which is identical to what we found from zero temperature limit of the thermo- 
dynamics (see (3.17)), the expression for the excitation energy turns into 

2JV-2 

E,^^ = e^i\r)- e,(Aj) . (4.14) 
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By using the integral equations for e^^ and one can show that 

e„(A) = . 

Therefore the excitation energy is exactly the same as for an excitation with only 
2N — 2 holes and no strings. The momentum for odd N is given by 

P = -— 2^ + -J- mod 27r , 

whereas for even N we find 

As we are only considering excitations with a finite number of particles in the ther- 
modynamic limit, the range of integers is microscopic [N — n <^ L) and the con- 
tribution of to the momentum goes to zero. Thus the momentum is identical to 
the momentum of an excitation with 2[N — 1) holes only. The spin quantum numbers 
are found to be 

= N-n = S . 

Excitations with an even number of holes and more than one string of length 
greater than 1 can be contracted in the same way as above. The result is that 
only the holes are dynamical, and the strings only affect the spin quantum number. 
This phenomenon is analogous to the situation for the Heisenberg magnet t^'*'"^^!. The 
interpretation of the spin- wave type excitations follows from the work of Faddeev and 
Takhtajan on the Heisenberg ferromagnet : in the spin sector of the m(2|2) model 
there exists one two elementary excitations of spin ib| and charge zero (which are 
called spinous), but they are "confined" and only their scattering states (triplet, 
singlet, ...) are observable. As an example we consider the sectors of 2 and 4 spinous 
respectively. For 2 spinous there is the state with N = 1 (of the first type of excitations 
considered above), which has spin S = 1 = and is the lowest weight state of the 
triplet. The other state is the one with N = 2 = n (of the second type of excitations 
considered above), which is the singlet state with S = = . Thus we obtain the 
correct irreducible representaions of the spin SU(2) for 2 particles. 

In the sector with 4 spinous we have the state with N = 2, which is the lowest 
weight state of the symmetric S = 2 representation, three states with = 3, n = 2 (as 
l-^al ^ 1), which lead to three = I triplets, one singlet with A^ = 3 = n, and another 
singlet with M-^ = N^^ + A^P - 4, Aff ^ = N^^ , = 2 (which is a states with 4 
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holes and 2 strings of length 2). Altogether these are 16 states in the representation 
|®|®|®| = 2®1®1®1®0®0 as expected. This completes our analysis for the 
Bethe states in the spin-wave sector. 

Each of the excitations constructed above is not only the lowest weight state of 
a spin SU(2) multiplet (as we took into account above), but of a m(2|2) multiplet, 
the dimension of which is given by (1.7) or (1.8). Thus we can construct further 
excitations with different quantum numbers by applying the m(2|2) raising operators. 
The energies of these excitations follow from the dispersion of the spin-waves and 
(1.9). We see that all excitations in the multiplet except the spin-SU(2) descendants 
acquire a gap with respect to the ground state for less than half filling. For example 
the state obtained by acting with Q_^ on the spin singlet excitation constructed 
above, carries quantum numbers = ^ and charge -\-e (the electron is taken to have 
charge — e) and has a gap of ^ -|- // > 0. 



Spectrum of Charge Excitations 

Excitations connected with charge rearrangements that do not involve local elec- 
tron pairs, are called holon-antiholon excitations and are completely analogous to 
the ones in the supersymmetric t-J model. Therefore we will restrict ourselves to a 
brief discussion. A holon-antiholon excitation corresponds to rearrangements in the 
A'^^^-distribution only. It is a particle-hole excitation, which is constructed by remov- 
ing one of the integers (corresponding to spectral parameter A^) from the ground 

state distribution and placing it outside the "Eermi boundary" — ^ (which is the 

maximal integer taken in the ground state). We denote the spectral parameter cor- 
responding to this newly "taken" integer by A^. Taking the thermodynamic limit we 
obtain the following set of coupled integral equations for the densities 

12 I 4 

,<..„,..) = _L £ cIA^^^jl^ ,,,A) - i*-,A." - A,, . 

Defining = —8{K — Kj^) and = ^ i_|_(a^-a)2 ? these equations can be written 
in matrix form (4.5), and following through the same steps as for the spin excitations 
above we find the excitation energy to be 
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/OO J-B TJ 

■OO J-B ^ 

(4.15) 

= E / £„(/^J V-^ +/^+i7 = 4^^(A,)-4^^(AJ>0. 

Here are the backllow distributions defined in (4.7). The contribution // + ^ is due 
to the lact that L J^^ dA p'^^^A) = N^^ — 1, as the particle with is outside the 
Fermi boundary. The excitation energy clearly breaks up into a contribution from the 
hole —e'--^\Af^) > and from the particle e'^-^^Ap) > 0. The particle carries no spin and 
charge +e, whereas the hole carries no spin and charge — e. Adding the particle at the 
Fermi boundary ±B turns the two-parameter family into a one-parameter one which 
would describe only a hole. At half filling B = the complete A*^^^ sea is empty and one 
can make neither particle-hole excitations nor an elementary particle excitationst^^^ 
with no spin (one gets scattering states of holons and spinous). Once again one can 
construct the complete u[2\2) multiplet from this excitation by acting with the raising 
operators. As the charge excitations all are spin singlets, the corresponding multiplet 
is of dimension 16[Nf^ -\- 1), where Nf^ -|- 1 is the dimension of the ?y-SU(2). However 
all these excitations have a gap for // < 0. For // = the localonic excitations become 
gapless (see below). 

Localonic Spectrum 

Let us now turn to the study of excitations involving local electron pairs. Such 
excitations are not present in the supersymmetric t-J model. We first turn to exci- 
tations involving A'^^^-strings. We consider the following situation 

Ml = iVf ^ + iVf ^, M„ = Vn > 2 , 
Mi'^ = N'f, M(2) = 1 , m(') = = Mf ) yq^p. 

This excitation corresponds to replacing p holes (empty sites) in the ground 
state with a A*^^^ string of length p. The excitation is spinless = and carries charge 
—2pe. It describes an excitation of p bound local electron pairs. The integers 
are distributed exactly like in the ground state, whereas the jump from half-odd 
integers to integers. The allowed range of integers for the A*^^^ string follows from 
(2.12) and is found to be \N^ \ < K^^^"^ ~ 2p). The spectral parameter describing the 
center of the p-string is determined by the integer via the equation 

27riV[ = ^ 0{^^) . (4.16) 

7=1 
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The integral equations for the A and A*^^^ distribution functions are found to be 



27r 1 + (A-A(i))2 ' 2ttL Vp' + (A(^) - k)^ 

If we define 



1 2p 



27rp2 + (A(i) - Kf 

the backfiow-distribution once again obey a matrix integral equation of type 
(4.8). The excitation energy follows from (2.13) to be 

E{k) = E{B) - E^siBo) = -2pfi + J2 dfi^ i^a 



1 r 2p 



Using the expression for the dressed energy of a A'^^^-string of length p given by 
(3.17) this can be written in the form 

E{k) = €^f{K) > -2iip > 0. 

We see that this kind of excitation has a gap for // < 0, which is the case we are 
currently investigating. This gap goes to zero as fi approaches zero. 

The dressed momentum of A'^^^-string excitations can be evaluated using (2.14) 
and (4.16) 



P{.) = f ( ±^ -nA= -J JA e{^) p[^\A) 



Note that the it sign is due to the fact that there exist two possibilities of 
distributing the integers = J^'^ it |. Let us concentrate on the case with the plus 
sign (the minus case can be treated in the same way). In the limit k, — )• oo, k,^ 
the momentum goes to zero as 



B 

d\ ■ 

K — A K 



P{.) I dA^ p[^\A) ~ ^Df' , 
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where Df^ = —j— is the density of holes in the ground state. We observe that zero 
momentum corresponds to infinite k. The excitation energy to leading order in ^ is 



E(k] 



Combining the asymtotic expressions for energy and momentum, we find a 
quadratic dispersion relation with a gap proportional to the magnitude of the chemical 
potential for small momenta 

E{P) = -2pfi + ^ + 0{P'), 
where Ai plays the role of a mass and in our approximation is given by 



M 



j%dK (A) 



> 



Here/_^ JA e(,')(A) is the energy of the particles in the A'^^^-ground state distribution. 



o 
I 



+ 

X 




Figure 3: Dispersion relations for a X'^^^ -string of length 1 
for ji = —0.1 and various values of U . 
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Figure 3 displays the dispersion relations for a localonic excitation with p = 1, 
i.e. for a A'^^^-string of length 1, for various values of U and fixed chemical potential 
/i = —0.1. The gap for this particular excitation is — 2/i, and we have displayed 
the whole range of momentum by taking in account both signs for the shift of the 
integers (the + ( — ) sign leads to positive (negative) momentum, of the same type 
of excitation). From Figure 3 we see that the allowed range of momentum tends 
to zero as /i + ^ approaches 21n(2), which corresponds to the ground state filled 
with single electrons (half-filled t-J ground state). On the other hand, as U goes to 
zero we observe an "extended Brillouin zone" - the single particle excitations appear 
multivalued in the interval [— 7r,7r], and only appear as single valued in an extended 
Brillouin zone scheme in a larger interval. This is a well know phenomenon, see e.g. 
[36]. Fxcitations of A*^^^ -strings of length 1 with zero momentum correspond to the 
action of the symmetry generator r]K As zero momentum corresponds to infinite 
spectral parameter A*^^^ there is no Pauli-principle for local electron pairs with zero 
momentum (which is in perfect agreement with the m(2|2) symmetry of the model). 

It is rather instructive to investigate the zero-density limit in some detail. This 
limit corresponds to taking // -|- ^ to zero, which is equivalent to i? — )• oo. In this 
limit it is possible to solve the coupled integral equations (3.18) analytically with the 
result e'--^^{A) = —j:^- Thus the excitation energy for a single localon (A'^^^-string of 
length 1) follows to be £'(A^^^) = e!|^^^(A^^^) = — 2/i -|- q_|_^^(2)-|2 • The state is described 
by the quantum numbers = 1 = Np Nj = 1^. Naively one might have expected 
to arrive at the state with one local pair over the empty vacuum. Although this is 
indeed an eigenstate of the hamiltonian, it is not given by the Bethe Ansatz in the 
BBFF representation, as was shown in [10]. Only lowest weight states of the m(2|2) 
algebra are Bethe Ansatz states, and the one-localon state is a m(2|2) descendant of 
the state with one spin down (in the BBFF representation). 



Electronic Spectrum 

The other type of excitation involving local electron pairs is the A*^^^ — A'^^^-string 
excitation. It will turn out however, that is does not describe a localonic excitation 
but rather carries the quantum numbers of an electron. It is constructed by choosing 

Ml = iVf ^ + iVf ^, M„ = Vn > 2 , 

^{1) ^ ^GS _ ^(1) ^ ^ ^ ^(2) = V? . 

Like for the ground state all vacancies for the integers I\ are filled symmetrically 

^In the FFBB representation (for which the structure of the string hypothesis is the same as for 
the BBFF representation, see below) it is given by the state with one A string of length 3, and 
one A*''-A*^' string. 
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around zero. The J^'s change from half-odd integer to integer. However, the number 
M^^"^ of integers also decreases by one as compared to the ground state, so that 
they are again distributed symmetrically around zero. The allowed range of integers 
for the A(^) - A^'^-string follows from (2.12) to be equal to lA^^I < |(iVf^ - 1). The 
rapidity k, describing the center of the A*^^^ — A*^^^ string is determined by its integer 
via 



27rK, = J2e{^)-J2 . (4.18) 



a=l 7=1 



Once again the backflow distributions obey a matrix integral equation of the 
form (4.8) where now the functions t/j are given by 

14 ,12 



27r4 + (A-/€)2 ' 27rl + (A-/€)2 

In terms of dressed energies the excitation energy is given by 



1 /.oo 4 1/"^ 2 

E(k) = U+— dXeJX) / dAe\'\A) -. 

SttI^ i^^4 + (A-/c)2 27rJ_s i ^ ^ + (A _ «)2 

Comparing this expression with the zero temperature limit of the TEA (3.17) 
we see that this is precisely equal to the dressed energy of a A'^^^-A'^^^ string, i.e. 



i?(«)^4^)(«) 



The dressed momentum can be determined using (2.14) and (4.18) 



/CO pB 
dx p^{x) e{^) - dA p['\A) e{.-A) . 
■oo J-B 

The quantum numbers carried by this type of excitation are S = ^ = and 
charge — e, so that it may be interpreted as an electron. The dispersion relation is 
plotted below in Figure 4 for a fixed chemical potential of // = — 0.1 and various values 
of U. Clearly there is always a gap, which tends to zero for // + ^ — )• 0. In this limit 
however the allowed range of momenta also goes to zero, so that the state cannot be 
excited. On the other hand, when // + ^ — )• 2 ln(2) (when the ground state is equal 
to the half filled ground state of the supersymmetric t-J model) the allowed range of 
momentum is [— 7r,7r]. 
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U=0.436 



Q I I I I I I I I I I I I I I I I I I I I I I I I 

-1 -.5 .5 1 

p 

' ex 

Figure 4: Dispersion relation for a A'^^^-A'^^) string. 



For small momenta the dispersion relation is quadratic and has a gap 
E{P) = U-eM-^+ ,, ,-51(0) , + 0{P') . 

Note that is the region under consideration U > (ei(0) < always), and the gap is 
positive. For momenta close to the Fermi momentum kp the dispersion is of the form 

E{P) ^U-{P- kpf + 0((P - kpf) . 

This shows that the gap approaches U at the Fermi boundary, which is in perfect 
agreement with Figure 4 (note that in Figure 4 we have subtracted — 2/i from the 
energy). 

This concludes our analysis in region IIL In summary we find the gapless excita- 
tion spectrum to be the one of the supersymmetric t-J model. However there exists 
a large number of new excitations with a gap. 
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4.2. Excitations in Region I 

We have shown in [10] that the ground state in region I, that is for zero chemical 
potential and attractive on-site interaction [/ < 0, is of the form (?y^)^|0), where 2N 
is the total number of electrons (note that we consider the ground state only for 
even total numbers of electrons, because it is otherwise trivially degenerate). The 
analysis of the excitation spectrum in this sector is most easily performed in the 
FFBB representation (see [10]). The reference state in this representation is the 
empty vacuum |0), which is much "closer" to the ground state than the reference 
state of the BBFF representation used above. The ground state is of a similar nature 
as the purely ferromagnetic ground state in the spin-| Heisenberg XXX model, as 
it is also an SU(2) descendant of the bare (empty) vaccum |0). Due to this nature 
of the ground state one would not expect the excitations to get a nontrivial dressing 
from the ground state electrons as is the case for antiferromagnetic systems. This is 
indeed the case as we will now proceed to show^°. In [10] we have constructed the 
Bethe Ansatz eigenstates of the hamiltonian, which are of the form 

m 

X{C.S\)\^)F^--K (4.19) 

8 = 1 

Here the C^,{\^) are creation operators with respect to the vacuum |0). A com- 
plete set of eigenstates can be obtained by acting with the m(2|2) raising operators 
on these states. We furthermore evaluated all commutators of the generators of the 
m(2|2) algebra with the creation operators C . As a special case we find 

h^c.(A)] = 0. 

This implies that 

ra ra 

iri^TX{CaS\MF''-''' = X{CaS\)(r,^mF^--^ . (4.20) 

8=1 8=1 

Recalling that a complete set of states is obtained by actiong with the m(2|2) rais- 
ing operators of the states (4.19) we arrive at the following expression for a complete 
set of excitations over the ground state {'f]^)^\0) 

m 

«(2|2) X ]lC\i\Wf\0)F^--^ , 

8 = 1 

Note that this would not be obvious if we operated in the BBFF representation. 
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where u[2\2) X denotes action by all raising operators of u[2\2) and rj. Now we 
see that due to (4.20) the only nontrivial part of the energies comes from the Bethe 
states, as action with the u[2\2) generators just shifts the energies according to (1.9) 
(and the hermitean conjugate equations). Let us now look at some excitations in more 
detail. It was shown in [10] that the Bethe equations for the BBFF and the FFBB 
representations are identical up to the substitutions Nf^ N^- and <^ N^. Thus the 
logarithmic form of the Bethe equations and the boundaries for the integers derived 
in section 2 for the BBFF representation also hold in the FFBB representation up to 
these substitutions. Fnergy and momentum in the FFBB representation are given by 



CO M<i') 

+ Mf)(/.-^ + /.) + M(^)(2/.-[/)-5] J]2/.p , (4.2i: 



where now 

n 

M(^) := M^l^ + 2M^^^ = iV^ + iV^ , 

M(2) := M^^^ + ^ ^Uf^ = . 

p 

An excitation of a bound state of n local electron pairs over the bare vacuum |0) 
(and by the above also over the true ground state (?7^)^|0)) is easy to construct as 
it only involves one spectral parameter of the zeroth level Bethe equations. Picking 
= 1 = Af, M^^"^ = = M2^\ and Afj^^ = 0, we find the following dispersion relation 
(for h = and // = 0) 

4 / P 
E(P) = -sin^ - 
^ ' n \2 

where P is the momentum of the excitation with respect to the ground state. The 
quantum numbers are = and charge —2n. By means of the supersymmetry we can 
construct electronic excitations as follows : consider an excitation with a single local 
pair over the ground state as constructed above (n = 1). Denote the wave function of 
this state by \loc). The state Q^\loc) has a dispersion relation E[P) = 4sin^ (-y) — ^, 
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spin —cr and charge — e, i.e. describes an electronic excitation over the ground state. 
Acting with on \loc) leads to a "three particle excitation" (one single electron and 
one local pair) with the same dispersion as Q^\loc). Of course there is an infinite 
number of other types of excitations, e.g. A'^^^-string excitations of any length. Our 
discussion above shows how all of these can be constructed in principle from the Bethe 
Ansatz and the u[2\2) supersymmetry. 



4.3. Excitations in Region II 

At /i = the compressibility becomes infinite, and the ground state density is no 
longer determined by the value of the chemical potential^'']. It is therfore necessary 
to work in the canonical ensemble and fix the density of particles D = . The 

ground state in region II (that is for /i = 0, < [/ < 41n(2)) is of the form {ri^'f\stJ), 
where \stJ) is the ground state of the supersymmetric t-J model at chemical potential 
/ijj = ^, and n is determined by the density imposed. This region is of particular 
interest, as the ground state in this region exhibits off-diagonal long-range order^''] 
and is thus superconducting. The ground state can also be determined from the zero 
temperature limit of the thermodynamic equations: for // = the dressed energies 
of localonic excitations become gapless, and due to the absence of a Pauli principle 
for local pairs at zero momentum there is an infinite number of states with the same 
energy as the ground state obtained by taking the limit // — )• from the left. Thus 
there is no unique ground state in the grand canonical ensemble, but by imposing a 
fixed density D we can single out a unique ground state, which is of the above form. 
The excitation spectrum now can be easily constructed by combining the results of 
sections 4.1. and 4.2.. By taking the limit // — )• of the spectrum determined in 
4.1. we obtain the excitations over the state |stJ), which by (4.20) readily translates 
into the excitation spectrum over the ground state {r]^)"'\stJ) . The result for < 
U < 41n(2) is as follows: the spin- wave and charge excitations are like for // < 
(gapless, linear dispersion around their Fermi surfaces). The localonic excitations 
are now gapless^ but still have a quadratic dispersion. The electronic excitation still 
has a gap proportional to U. The quadratic behaviour of the gapless dispersion of 
the localons clearly indicates non-conformal properties of the model. The methods 
of conformal field theory^"^''"'*^] can therefore not be applied directly to evaluate the 
critical exponents. However, by using the algebra of the m(2|2) generators we have 
succeeded to express correlators in the m(2|2) model is terms of correlators in the 
supersymmetric t-J model. This enables us to evaluate the critical exponents using 
the conformal field theory results for the t-.J model^'*"^''*'*]. 

Like for // < every Bethe Ansatz excitation provides the lowest weight state of 
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a whole u[2\2) multiplet, the structure of which is given by (1.7) and (1.9). 
For U > 41n(2) \stJ) is the half-filled ground state of the supersymmetric t-J model, 
for which the excitation spectrum is rather different as compared to less than half- 
fiUingt^^l. We leave the discussion of the region for a future publication. 



4.3. Excitations in Region V 

Excitations in region V can be constructed the same way as in region I. The 
main difference is that the chemical potential in region V is not zero, and that only 
the action of the u[2\2) raising operators has to be taken into account. 
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